Ground-state properties of one-dimensional matter and quantum dissociation of a 

Luttinger liquid 
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Motivated by emerging experimental possibilities to confine atoms and molecules in quasi-one- 
dimensional geometries, we analyze ground-state properties of strictly one-dimensional molecular 
matter comprised of identical particles of mass m. Such a class of systems can be described by an 
additive two-body potential whose functional form is common to all substances which only differ 
in the energy e and range / scales of the potential. With this choice De Boer's quantum theorem 
of corresponding states holds thus implying that ground-state properties expressed in appropriate 
reduced form are only determined by the dimensionless parameter Aq ~ h 2 /ml 2 e measuring the 
strength of zero-point motion in the system. The presence of a minimum in the two-body interaction 
potential leads to a many-body bound state which is a Luttinger liquid stable for a not very large 
Ao. As Ao increases, the asymmetry of the two-body potential causes quantum expansion, softening, 
and eventual evaporation of the Luttinger liquid into a gas phase. Selecting the pair interaction 
potential in the Morse form we analytically compute the properties of the Luttinger liquid and 
its range of existence. We find that as Ao increases, the system first undergoes a discontinuous 
evaporation transition into a diatomic gas followed by a continuous dissociation transition into a 
monoatomic gas. Two-body potentials of molecular systems can be successfully fitted into the Morse 
form thus allowing determination of the quantum parameter Ao and the state of matter of substance 
in question. In particular we find that spin-polarized isotopes of hydrogen and 3 He are monoatomic 
gases, 4 He is a diatomic gas, while molecular hydrogen and heavier substances are Luttinger liquids. 
We also investigate the effect of finite pressure on the properties of the liquid and monoatomic gas 
phases. In particular we estimate a pressure at which molecular hydrogen undergoes an inverse 
Peierls transition into a metallic state which is a one-dimensional analog of the transition predicted 
by Wigner and Huntington in 1935. 

PACS numbers: 68.65.-k, 61.46. +w, 71.10.Pm, 05.30.-d 



I. INTRODUCTION 

Prediction of the ground-state properties of a con- 
densed many-body system of identical particles starting 
from microscopic two-body interactions is only meaning- 
ful if the outer electronic shells of underlying atoms or 
molecules in bound phase are not very different from their 
free state counterparts [1]. If the two-particle poten- 
tial has a functional form common to a family of sub- 
stances (for example, of the Lennard- Jones type), then 
the properties of all the members of the family can be 
related. This conclusion pioneered by De Boer and col- 
laborators, commonly referred as the quantum theorem of 
corresponding states, was originally applied to predict the 
properties of 3 He [2] before it had become experimentally 
available. Later Anderson and Palmer [3] and Clark and 
Chao [4] used the same approach to estimate the proper- 
ties of zero-temperature nuclear and neutron-star matter 
from those of laboratory substances. 

The goal of this work is to conduct a similar program 
in a strictly one-dimensional case. There are several rea- 
sons, both of fundamental and practical nature, why it 
is important to understand this problem. 

First, it is well-known that for ordinary substances 
zero-point motion is of crucial importance only for the 
lightest elements such as helium isotopes as well as spin- 
polarized isotopes of hydrogen. The ground state of all 



heavier elements (and molecular hydrogen) is crystalline 
and to large extent classical. The one-dimensional case 
is qualitatively different: regardless of the particle mass 
zero-point motion destroys the long-range crystalline or- 
der - a situation that is closely analogous to the destruc- 
tion of the long-range order by thermal fluctuations in 
classical two-dimensional systems of continuous symme- 
try [5]. As a result, the only possible many-body bound 
state in one dimension is a harmonic or Luttinger liquid - 
a uniform density condensed phase with algebraically de- 
caying density correlations [6] . This decay, characterized 
by a nonuniversal exponent, is slower than the exponen- 
tial fall-off of density correlations in conventional fluids. 
As the degree of zero-point motion increases, the many- 
body bound state can disappear through a transition that 
has no analog in the three-dimensional world: since the 
Luttinger liquid phase is more correlated than conven- 
tional fluids and less correlated than standard crystals its 
dissociation combines qualitative features of both labo- 
ratory melting and evaporation at once. 

Second, in one dimension the difference between 
fcrmions and bosons is not very significant as we cannot 
go from one configuration to another with exchanged par- 
ticles without bringing the particles in contact at some 
intermediate step. Then short-distance repulsion among 
bosons has the same effect on density correlations as the 
Pauli principle for fermions. On the other hand, zero- 
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temperature properties of three-dimensional matter with 
an interaction pair potential of the Lennard-Jones form 
are sensitive to the statistics of the underlying particles 
[7]-_ 

Since experimental discovery of carbon nanotubes in 
1991 [8] studying the properties of one-dimensional sys- 
tems became especially important. In addition to their 
unique transport, mechanical and chemical properties 
[9] , bundles of carbon nanotubes can play a role of one- 
dimensional hosts for foreign atoms that can find them- 
selves bound in the interstitial channels or inside the 
tubes [10]. One of the interesting potential applications 
of these systems includes storage devices for molecular 
hydrogen in fuel cells [11]. 

Recently the quasi-one-dimensional regime has been 
also realized for Bose-condcnsatcs of alkali vapors both 
for repulsive [12] and attractive interactions [13]. These 
systems which are relevant for atom interferometry [14] 
have an additional flexibility as the strength and sign of 
two-body interactions can be magnetically tuned. 

In both of these experimental examples the basic model 
is a zero-temperature one-dimensional many-body sys- 
tem of fermions or bosons with pairwise interactions. 
There were several attempts in the past to study this 
problem: 

(i) Diffusion Monte Carlo studies predicted that at 
zero temperature both one-dimensional 4 He [15,16] and 
molecular hydrogen [17] form weakly-bound liquids. 

(ii) These conclusions were supported by variational 
studies based on the Jastrow-Feenberg wave function [18] 
where additionally it was argued that in one dimension 
the many-body bound state exists only for those systems 
which have a dimer, i. e. a two-body bound state. One 
of the findings common to Refs. [16-18] is the prediction 
of a high-density liquid-solid phase transition in which 
a standing density wave sets in. We note however that 
the existence of such a one-dimensional solid contradicts 
the quantum version of the Mermin-Wagner-Hohcnbcrg 
theorem [5]. 

(iii) A direct variational treatment based on a Gaussian 
wave function was performed in Ref. [19] where it was 
assumed that the particles form a one-dimensional chain 
with the Morse potential interaction [20] between near- 
est neighbors. Although it was shown that the chain re- 
mained stable for not very strong quantum fluctuations, 
the accuracy of the method, the nature of the condensed 
phase, the role of dimerization, and the implications for 
real systems were not addressed. 

In this paper the problem of the ground-state prop- 
erties of a one-dimensional many-body system is re- 
examined for the case when the two-body interparticle 
interaction can be approximated by the Morse potential 
[20]. Since it involves three parameters, the Morse po- 
tential is more flexible in describing real systems as com- 
pared to the two-parameter Lennard-Jones potential. At 
the same time the quantum theorem of the corresponding 



states [2] still holds in this case. In addition the prob- 
lem of the Morse dimer is exactly solvable [20] ; below we 
will also show that analytical progress is possible in the 
many-body case, and the accuracy of our results can be 
assessed. 
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FIG. 1. The Silvera-Goldman interaction potential be- 
tween two hydrogen molecules (crosses) and its approxima- 
tion by the Morse potential (solid curve) . To demonstrate the 
strength of the overlap repulsion, logarithmic scale is used on 
the positive energy axis. The details of the fit are explained 
in the text. 

The organization of the paper is as follows. In Section 
II we set up the problem in general terms with the quan- 
tum theorem of the corresponding states [2] as a guide 
for possible outcomes. The parameters of the Luttingcr 
liquid arc computed in Sections III-V. In Section III wc 
describe general properties of the Luttinger liquid and 
outline the main idea of the calculation. Zero-pressure 
analysis is carried out in Section IV. The main tool here is 
a combination of variational and renormalization-group 
treatments. As a by-product we also solve the problem of 
quantum Brownian motion in the Morse potential, and 
show that it exhibits a localization-delocalization tran- 
sition. Finite-pressure variational analysis is conducted 
in Section V. Section VI is dedicated to the discussion 
and applications of our results to various molecular sub- 
stances. In particular we estimate a pressure at which 
one-dimensional molecular hydrogen undergoes a transi- 
tion into a metallic state. 

II. FORMULATION OF THE PROBLEM 

Our starting point is the Euclidian action for N iden- 
tical fermions or bosons of mass m with pairwise inter- 
actions 

r m N At n 

i=l i<j 
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where r is the imaginary time variable, Xi are particle 
positions, and V(h) is the pair interaction potential. The 
properties of the function V(h) can be summarized as 
follows. At large separation h interparticle interaction 
is dominated by weak rapidly decaying Van der Waals 
attraction, while at short distances there is a very strong 
overlap repulsion [1]. As a result, the pair potential V(h) 
has an asymmetric minimum at some intermediate h. As 
a typical example Fig. 1 shows the semi-empirical Silvera- 
Goldman potential [21] between two H2 molecules which 
is extensively used in computations of the properties of 
molecular hydrogen. 

Assume that the pair potential has the form 



V(h) = eU[(h-H )/l], 



(2) 



where e is the energy scale of the potential, I is the po- 
tential range, H is a length scale, and U (y) is a function 
common to a family of substances. Introducing dimen- 
sionless position and time variables, = Xi/l, t — er/K, 
transforms the reduced action S/% into 



S 

h 



2—1 i<j 



where Qo = Hq/1. The quantum theorem of the corre- 
sponding states [2] then directly follows from representa- 
tion (3): the energy per particle E* measured in units of 
e is only determined by the dimensionless parameter 



A 



nl^me) 1 / 2 '' 



(4) 



the form of the function U in (2), and the particle statis- 
tics: 



E* = E*(Xq, statistics) 



(5) 



Similar statements can be made about the reduced rel- 
ative equilibrium length per particle Q — Qo (Q is the 
one-dimensional version of volume per particle measured 
in units of I), and other quantities of interest. Apart from 
numerical factors (introduced for convenience), the quan- 
tum parameter A (4) is identical to De Boer's number 
[2]: its square is proportional to the ratio of the zero- 
point energy of particle of mass m localized within range 
I to the typical potential energy e. Therefore as A in- 
creases away from its classical value Ao = 0, the strength 
of zero-point motion increases. 

In the presence of several competing phases the func- 
tion E* in (5) (and other properties) will have sev- 
eral branches; the branch with lowest E* singles out 
the ground-state of the system. When two different 
branches cross, the ground-state changes via a first-order 
phase transition. Each branch of E*(Xq) is an analyti- 
cal function of its argument except possibly at isolated 
points where critical phenomena take place. One obvi- 
ous branch of E* corresponds to a monoatomic gas which 



must become the ground state of the system at suffi- 
ciently large Ao- Then all the particles are infinitely far 
apart from each other, and thus -Bm 0no (Ao) = which we 
select to be the reference point for the energy. 

In what follows we select the pair interaction potential 
in the Morse form [20]: 



V(h) = -Ae- h ' 1 + Be- 2h ' 1 

= e[-2e~ {h ~ Ho)/l + e - 2 (>*- H o)/i^ 



(6) 



where A and B arc the amplitudes of the attractive 
and repulsive parts of the potential, respectively, Hq = 
lln(2B/A) is the location of the minimum of (6), while 
e = V(H ) = A 2 /AB is the depth of the potential well. 
The second representation of (6) shows explicitly that 
the Morse potential conforms to the general form (2). It 
is physically reasonable to require that the zero of (6) is 
located at positive h which implies B > A. 

Similar to the applications of the Lennard-Jones po- 
tential to laboratory molecular systems [1], Eq.(6) should 
not be taken too literally as really describing two-particle 
interactions. The only reason behind our choice (6) is the 
possibility of analytic progress. 

For two particles interacting according to (6) the 
ground-state energy is exactly known to be E2 = — e[l — 
h/2l(me) 1 / 2 } 2 [20]. This implies that in the many-body 
case one of the possible phases of the system is a diatomic 
gas (a collection of infinitely far separated dimers) with 
the reduced energy function 



E dimer(^o) - ~ ""/f^ 



(7) 



valid for Ao < A02 = V2/tt; the factor of 1/2 accounts 
for two particles in the dimer. As Ao approaches A02 
from below, the dimer size diverges, and at Ao = A02 a 
second-order dissociation transition into the monoatomic 
gas discussed earlier takes place. The asymmetry of the 
interaction potential is responsible for the disappearance 
of the two-body bound state for sufficiently strong zero- 
point motion. 

The diatomic gas might be the ground-state of the sys- 
tem for intermediate Ao but for sufficiently small Ao a 
condensed phase must have the lowest energy. For molec- 
ular systems in general [1] pair interactions decay rapidly 
with interparticle separation. As a result, the physics of 
the condensed phase is dominated by nearest-neighbor 
interactions. Therefore in what follows in describing the 
one-dimensional condensed phase we restrict ourselves to 
nearest-neighbor interactions. Corrections coming from 
ignoring distant neighbors will be marginally small pro- 
vided the interaction range I is substantially smaller than 
average interparticle spacing (bond length). 

In the classical limit, A = 0, the ground-state of the 
system is a one-dimensional crystal of particles with lat- 
tice spacing H given by the minimum of the two-body 
Morse potential (6). Indeed, the energy per particle for 
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the crystal, — e, is twice as negative as that for the di- 
atomic gas. 

A condensed phase is also the ground-state of the sys- 
tem for a not very large Ao- The quantitative theory of 
the properties of this phase is developed below. 



III. LUTTINGER LIQUID PHASE AND ITS 
PROPERTIES 



When zero-point motion is present, the asymmetry of 
the pair interaction potential about its minimum causes 
quantum expansion. As a result the average bond length 
H in the condensed state exceeds its classical, Ao = 0, 
counterpart Hq. The low-energy dynamics of the system 
is described by the harmonic action [6] 



Sha 



(8) 



where u(x, r) is the particle displacement field, p = m/H 
is the mass density, and c is the sound velocity [1] 



c2 _ H 2 _d 2 E{h = H) 



m 



8h 2 



(9) 



where E{h) is the ground-state energy per particle as a 
function of (one-dimensional) volume per particle h, and 
the derivative is evaluated at the equilibrium interparticle 
spacing H. The function E(h) can be also viewed as 
an effective pair interaction renormalized by zero-point 
motion away from its classical form (6). 

The Feynman path integral formulation of quantum 
mechanics [22] allows us to view the action (8) as a 
Hamiltonian for a classical two-dimensional crystal of 
line objects (world lines of underlying particles) running 
in the imaginary time direction. In this correspondence 
zero-point motion plays a role of thermal fluctuations. 
But this is exactly the context of applicability of the 
Mermin-Wagncr-Hohcnberg theorem [5]: if n(x, r) is the 
instantaneous number density, then long-wavelength low- 
energy quantum fluctuations captured by the action (8) 
destroy long-range positional order of the particles. The 
only allowed many-body bound state is a uniform density 
phase, < n(x,r) >= H~ x (<> stands for the expecta- 
tion value), with algebraic decay of density correlations 
[6]: 



< n(x, r)n(0, 0) > -H~ 2 oc 
where the exponent g is given by 

9 = ^h 2 



cos(2irx / H) 

(X 2 + C 2 T 2 )9 ' 



(10) 



(11) 



The large distance/time behavior of the density-density 
correlation function (10) is the hallmark of the Luttingcr 



liquid. In order to compute the range of existence of the 
Luttinger liquid phase, its bond length H, sound veloc- 
ity c, and correlation exponent g we need to go beyond 
harmonic approximation. 

Our calculation relies on the assumption underlying 
the harmonic description (8) (and to be verified later) 
that as long as the Luttinger liquid is stable, the ra- 
tio of the typical fluctuation of the bond length to the 
bond length itself remains small despite individually both 
these quantities are increasing functions of Ao (4). Then 
every bond of the system can be viewed as a quantum- 
mechanical degree of freedom subject to the external po- 
tential V(h) and placed in contact with a bath of har- 
monic oscillators (8) corresponding to the rest of the sys- 
tem. The single bond dynamics is thus given by the ac- 
tion: 



Sbo 



J dxdr[{ 



du^ 9 , du 
+C{ d~x 



) 2 } + J drV(h) (12) 



bath 



The first integral is over all positions and times except 
for a small region near x = where the bond in question 
is located. This separates the system into two pieces, so 
that at all times the displacement field u is discontinuous: 
u(x — +0,t) — u(x = — 0,t) = h(r) — H - The coupling 
between the segments of the system joined at the bond is 
given by the foil pair potential V(h), i. e. it goes beyond 
harmonic approximation. 

The action (12) has been previously introduced in Rcf. 
[23] to describe tunncling-assisted fracture of a stretched 
one-dimensional chain. More generally this type of ac- 
tion corresponds to the Caldeira-Leggett model of cou- 
pling between a quantum-mechanical degree of freedom 
and an environment modeled by a reservoir of harmonic 
oscillators [24]. Since the last interaction term in (12) is 
restricted to a single spatial point, the bath degrees of 
freedom can be integrated out away from the bond with 
the result [23]: 



S = 



pc 



cLu 
27 



drVQi), 



(13) 



where the subscript is dropped for brevity and the 
Fourier transform of the bond length field h{uo) = 
J h(r) exp(— wr)dr has been introduced [25]. The fre- 
quency cutoff ojd setting the limits of integration in the 
first kinetic term of (13) is the one-dimensional Debye 
frequency ojd = nc/H - the vibrational spectrum of the 
system is approximated by the Debye model. 

The unusual \u>\ dependence of the kinetic term of the 
action (13) is due to the many-body nature of the bond 
dynamics, and can be understood heuristically by notic- 
ing that if the bond length oscillates with frequency u, 
then during one oscillation period 27r/|w| this disturbance 
propagates in both directions away from the bond a dis- 
tance of order c/|w|. Therefore the standard kinetic en- 
ergy density, proportional to pu) 2 should be multiplied by 
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the size of the region c/\u>\ affected by the motion thus 
reproducing the pc\u>\ term of (13). 

The action (13) allows us in principle to compute how 
the bath degrees of freedom renormalize the properties 
of a given bond. The nonanalytic \u>\ dependence in (13) 
guarantees that the bond cannot renormalize the bath 
oscillators (whose properties are accumulated in the pc 
combination). This observation combined with the fact 
that all the bounds of the Luttinger liquid are equivalent 
provides us with a prescription on how to use (13) to 
solve the problem we are interested in: 

The parameters p and c of (13) should be selected as 
initially unknown but fully renormalizcd properties of the 
Luttinger liquid. The reservoir degrees of freedom will 
renormalize the microscopic pair interaction V(h) into a 
form which we will require to be identical (in the har- 
monic limit) to the rest of the chain. This will determine 
the parameters of the Luttinger liquid and guarantee that 
the treatment is insensitive to the choice of the bond. 

Since the action (13) describes the dynamics of an ar- 
bitrary single bond of the system, it can be directly used 
to compute the energy per particle of the original many- 
body problem. 



IV. ZERO-PRESSURE ANALYSIS 

For a quantitative analysis we use Feynman's [26] vari- 
ational principle for the ground-state energy: 



E < Ei = E + (T/h) <S-S > 



(14) 



where T is the temperature, and %/T has a meaning of 
the system size in the r direction; the T = limit will be 
taken at the end. The notation <> denotes an expecta- 
tion value computed using an arbitrary reference action 
So, and Eq is the ground-state energy corresponding to 
So- 

This method has been remarkably successful in ana- 
lyzing the roughening phase transition [27], multilayer 
adsorption phenomena [28], wetting transitions [29], the 
problem of quantum Brownian motion in a periodic po- 
tential [30], the Coulomb blockade problem [31], and a 
variety of problems of quantum mechanics and quantum 
field theory [32]. 

It is physically reasonable to select the trial action So 
in a Gaussian form similar to that in [28], [29], and [31]: 



So 



pc 
4 



2^ 



co\\h(u)\ 2 + 



K 



I 



dr(h-H) 2 }, (15) 



where two variational parameters which include the fa- 
miliar bond length H and a new parameter K (control- 
ling the extent of fluctuations about H) are selected to 
minimize E\ in (14). The stiffness parameter K has 



a meaning of the curvature of the effective pair poten- 
tial evaluated at its minimum H, and appearing in (9), 
K = d 2 h E(h = H). 

Introducing / = h — H, the deviation of the bond 
length away from its equilibrium value H, the reduced 
root-mean-square (rms) fluctuation, /* = (< f 2 > 
/I 2 ) 1 / 2 , can be computed with the help of (15) as 



r = (2A) i /2i n v2 (1+7 -i )i 



where 



A 



■npcl 2 



(16) 



(17) 



is the dimensionless parameter quantifying the strength 
of zero-point motion in the Luttinger liquid, and 



7 = 



2K 

OJDpC 



(18) 



is the dimensionless counterpart of K. For the Morse 
pair interaction (6) the classical sound velocity is Co = 
(H /l)(2e/m) 1 ^ 2 . It is then straightforward to verify that 
in the classical limit, h — > 0, the quantum parameter A 
(17) reduces to De Boer's number Ao (4). 

In the Luttinger liquid the strength of zero-point mo- 
tion is characterized by the correlation exponent g (11) 
which is related to A (17) by 



Xn 2 /Q 2 



(19) 



where Q = H/l is the reduced bond length. This re- 
lationship demonstrates that if the interaction range I 
is known, then measuring the density-density correlation 
function (10) will allow us to compute A, and thus the 
rms fluctuation of the bond length (16). 

Using (15), (16), and (6) the reduced upper bound 
E* = E\/e entering (14) can be computed as 



£T(7, Q) = v- 1 ln(l + 7 ) - 2e Q °- Q (l + 7" 1 ) 



1\A 



where 



+ e 2(Qo-Q) (1+7 -l ) 4A j 



2ne 

Thud 



(20) 



(21) 



is related to the reduced Debye temperature as 9* = 
hujrj/t = 2n/v. 

The expression for E* should be minimized with re- 
spect to 7 and Q, and in case of multiple solutions the 
one minimizing (20) must be selected. 

Minimizing E* with respect to Q we arrive at the ex- 
pression for the reduced bond length which accounts for 
quantum expansion 



Q = Qo + 3Aln(l+ 7 - 1 ) 



(22) 
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Substituting this back into (20) the expression for the 
reduced energy E* can be written as 



(23) 
(24) 



E*( 1 ) = v- 1 \n(l+ 1 )-(l + 1 - 1 )- 2 \ 
Minimizing (23) with respect to 7 we find 

7 = 2v\(l + 7~ x r 2A 

A. Approximate solution: A = Ao 



First we look at a simplified version of the original 
problem when only one bond of the system is subject 
to the Morse potential (6) while the rest of the chain is 
harmonic. This situation is described by the action 



S = 



poco f duo 



4 J 2tt 



+ J d T {-Ae- h ' 1 +Be- 2h ' 1 ), 



(25) 



where po = tu/Hq, cq — (HQ/l)(2e/m) 1 ^ 2 , and ojdo = 
ttcq/Hq assume their classical values. The parameters A 
and B are selected so that in the harmonic approximation 
the Morse bond is identical to the rest of the chain. 

If the imaginary time coordinate r is viewed as a fic- 
titious space variable then (25) can be recognized as an 
effective Hamiltonian defining a classical statistical me- 
chanics problem [22]. If the bond field h is identified with 
an interface height and zero-point motion with thermal 
fluctuations, then this problem is a one-dimensional ana- 
log of the critical wetting problem [29], [33]. 

For finite Ao zero-point motion softens and lengthens 
the Morse bond while leaving the harmonic part of the 
chain intact. This is described by Eqs.(22)-(24) with A = 
Ao and v = v = 1/(tt\ ). Specifically, Eqs.(23) and (24) 
turn into: 



£T(7)=7r[Aoln(l+7)- J] 



7= -(1 + 7- 1 ) 

7T 



-2A 



(26) 



(27) 



First we note that 7 = is always a solution to (27) 
with E* = which describes two segments of the chain 
infinitely far apart from each other. As Ao increases away 
from the classical value Ao = 0, the parameter 7 in (27) 
monotonically decreases from 7 = 2/t: vanishing at Ao = 
1/2. For Ao > 1/2 only 7 = solves Eq.(27). When A 
approaches the critical value of 1/2 from below we have 



7 ■ 



ln(x/2) 
e 1-2*0 



(28) 



Correspondingly, the reduced energy (26) monotonically 
increases with Aq: in the classical limit, Aq — > 0, it rises 



linearly with Ao accordingto E* = — 1+\q[tt ln(l+2/7r) + 
2 ln(l + tt/2)] ~ -1 + 3.4361A . Since A ~ l/m 1 / 2 , Eq. 
(4), the energy E* is not analytic as 1/m — > 0, which 
is the relevant quantity [3]. This is expected because 
the crystal (Ao = 0) is qualitatively different from the 
Luttinger liquid (finite A ). Upon approaching A = 1/2 
from below the reduced energy vanishes as 



E* 



In(^/2) 

2A )e ^ 



(29) 



As A approaches 1/2 the bond length diverges and for 
Ao > 1/2 the two segments of the chain are infinitely far 
apart from each other. 

The behavior of the reduced bond length and its rms 
fluctuation just below the phase transition can be found 
by combining (16) and (22) with (28): 



Q 



31n(7r/2) 
2(1 -2A )' 



\n 1/2 (n/2) 
(I-2X0) 1 / 2 



(30) 



We note that although both of these quantities diverge 
upon approaching the phase transition, the relative fluc- 
tuation, f*/Q, vanishes. Thus fluctuating segments of 
the chain never overlap and our treatment is consistent. 

The essential singularities (28) and (29) at A = 1/2 
as well as the divergences (30) parallel those found in the 
context of wetting transitions [29]. 

The most valuable feature of the variational approach 
is its nonperturbative nature. The accuracy of varia- 
tional predictions depends on how close is the variational 
guess to the physical reality. For the problem defined by 
the action (25) the accuracy of our approach can be as- 
sessed and the special role played by Ao = 1/2 can be 
re-established by using a renormalization-group method. 



B. Perturbative renormalization-group treatment: 

A = Ao 

Following the argument of Brczin, Halpcrin, and 
Leibler originally given in the classical context of wet- 
ting transitions [33] we treat the Morse potential term 
in (25) as a perturbation. Then the lowest-order 
renormalization-group equations have the form: 



din ,4 
dln(uj D aC) 



(Ao + 1), 



dlnB 



(4A + 1), (31) 



where £ is the running scale in the r direction, and the 
equations describe how the Morse parameters A and B 
renormalize upon (i) successive integration out of high- 
frequency modes (first terms) followed by (ii) scaling 
transformation which restores the cutoff to its original 
value (second terms). Instead of following separate evo- 
lution of the coefficients A and B, it is more appropriate 
to look at the depth of the Morse well (6) e = A 2 /4B 
which is also proportional to the potential curvature at 
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its minimum. For its dimcnsionless counterpart v (21) 
Eqs.(31) imply 

dv 



= (1 - 2X )v 



(32) 



rfln(a)DoC) 

The flow diagram corresponding to (32) is sketched in 
Fig. 2 where we also show the locus of initial conditions 
v = I/ttAq of the model (25). 



2/n — 



A A A 
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FIG. 2. The flow diagram of the model (25) . The arrows in- 
dicate the direction of the flow and the dashed line vo = l/ir\o 
is the locus of initial conditions. The stable part of the v = 0, 
Ao > 1/2 fixed line corresponds to two segments of the chain 
infinitely far away from each other. For A < 1/2 the pa- 
rameter v is relevant and the Morse bond joins the segments 
together. 

First we note that verticality of the flow lines is a rigor- 
ous property of the model (25) as the poco combination 
(and thus Ao) does not renormalize to any order in v. 
There are clearly two regimes separated by a phase tran- 
sition at the critical value Ao = 1/2. For Ao > 1/2 the 
parameter v flows to zero which means that the chain is 
broken into two infinitely separated segments. On the 
other hand, for Ao < 1/2 the parameter v grows un- 
der renormalization eventually leaving the perturbative 
regime v <C 1. This means that the Morse bond holds 
both segments of the chain together. Upon approach- 
ing the critical value Ao = 1/2 from below there is a 
divergent time scale £ (analog of a correlation length in 
standard critical phenomena) which can be found from 
the condition i>(£) ~ 1: 



J D0 



(33) 



The one-dimensional nature of the problem then im- 
plies that the reduced curvature 7 (18) of the renor- 



malized potential vanishes as £ . Similarly the criti- 
cal behavior of the reduced bond length Q can be found 
as Q ~ ]n[B(£)/A(£)] ~ (3/2) Ln(w D0 £)> while its rms 
fluctuation is /* ~ ln 1/2 (w D0 0- With £ given by (33) 
and logarithmic accuracy these results coincide with their 
variational counterparts, Eq.(28) and (30). If we set the 
bare parameter vq in (33) at 2/-7T (the crossing of the locus 
of initial conditions v = l/7rAo and the A = 1/2 line) 
the renormalization-group results would become identical 
to those of the variational approach. 

From perturbative renormalization-group treatment 
alone we would not be able to make reliable statements 
about the phase transition at Ao = 1/2 as the analysis 
is valid for wq < 1 while for the problem in question 
(25) one has v n = 2/n ~ 0.64 < 1 - it is on the border 
of applicability of perturbative theory. However com- 
bining the above results with the nonperturbatve vari- 
ational analysis makes a strong case. Since the latter 
produces the same answers in the region of parameters 
where renormalization-group results are less certain, we 
argue that the variational solution of (25) is very accurate 
in the range of Ao between the classical limit Ao = and 
the dissociation transition Ao = 1/2 which is described 
exactly. 

As a side observation we note that the problem (25) is 
most likely to be related to that of quantum Brownian 
motion of a particle in a periodic potential [30] . Superfi- 
cially the only similarity between the two is that (25) can 
be also viewed as describing quantum Brownian motion 
in the Morse potential which has no periodicity. 

The similarity between the problems becomes notice- 
able if one inspects their treatments. A comparison shows 
that our expression for variational energy (23) is identi- 
cal to that of the periodic version of the problem with 
the amplitude of the periodic potential proportional to 
our parameter v (21) and our A = A corresponding to 
l/2ct of Fisher and Zwerger [30]. With this identifica- 
tion renormalization-group equation (32) coincides with 
its periodic counterpart [30]. Both problems have der- 
ealization transitions of the same universality class driven 
by zero-point motion. 

We remind the reader that the action (25) is an ap- 
proximation to the original problem of the ground-state 
properties of the Luttinger liquid - only one bond is sub- 
ject to the Morse potential while the rest of the system 
is purely harmonic. This chain is obviously stiff er than 
the original system - a smaller level of zero-point motion 
will be necessary to cause dissociation of the Luttinger 
liquid. Therefore the results derived for the model (25) 
imply that zero-pressure Luttinger liquid phase cannot 
exist for A > 1/2. 

We also note that the critical value Ao = 1/2 is larger 
than the dimer dissociation threshold A02 = V2/n — 
0.4502. This can be understood qualitatively by notic- 
ing that the dynamics of the Morse bond joining two 
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half-infinite harmonic segments is more inertial (and thus 
more classical) than that of the Morse dimcr. Therefore 
with the same underlying particles a weaker level of quan- 
tum fluctuations (i. e. smaller A ) suffices to break the 
dimcr. 



C. Accurate solution 

In order to compute the properties of the Luttingcr 
liquid more accurately we have to impose the condition 
that all the bonds of the chain are equivalent. Then using 
definitions of De Boer's number Ao (4) and its Luttingcr 
liquid counterpart A (17) the parameter v (21) can be 
calculated as v = X/ttXq; for A = A it reduces to v = 
1/ttXq previously used in approximate treatment of the 
problem. Similarly the reduced curvature of the effective 
pair potential (18) can be computed with the help of 
Eq.(9) with the conclusion that 7 = 2/n. This is exactly 
what was previously found in the approximate analysis 
in the classical limit, A = 0, when indeed all the bonds 
of the chain are equivalent. 

Substituting v = X/irX 2 ,, and 7 = 2/n in Eq.(20) we 
arrive at the expression for the reduced energy as a func- 
tion of dimensionless "volume" per particle Q: 



E*(Q) = (ttA 2 /A) ln(l + 2/tt) - 2e Q °- Q+xla{1+7x/2) 



+ e 



2(Qo-Q)+4A ln(l+ir/2) 



Similarly Eqs.(16), (22)-(24) transform into 
/* = (2A) 1 / 2 ln 1/2 (l + 7 r/2) 

Q = Qo + 3Aln(l + 7r/2) 



E* LL 



-2Aln(l+V2)[ 7rAm(1 + 2/7T)-l] 

Ao = Ae- Aln ( 1+ "/ 2 ) 



(34) 



(35) 



(36) 



(37) 



(38) 



Eqs.(34)-(37) are the main results of this Section. We 
note that in view of Eqs.(10) and (19) Eqs.(35)-(37) give 
the dependence of the reduced rms bond length fluctu- 
ation /*, bond length Q, and energy per particle of the 
Luttinger liquid E* LL on the experimentally measurable 
quantum parameter A. 

The dependence of De Boer's number Ao on A is given 
by Eq.(38); we arc interested in the inverse dependence, 
A(Ao). 

The right-hand side of (38) has a maximum at A s = 
Ln _1 (l + 7r/2) ~ 1.0591 of magnitude A 0s = e- 1 ln~ 1 (l + 
7t/2) ~ 0.3896. For A < A 0s Eq.(38) has two roots for 
A but only the smaller one is physical. For A = X 0s 
these two roots coincide, and for A > A 0s Eq.(38) has 
no solutions - Luttinger liquid phase is no longer stable. 
We note that the Luttinger liquid cannot sustain the level 



of zero-point motion stronger than that corresponding to 
Xq s ~ 0.3896 which is smaller than the dimer dissociation 
threshold A 02 = V2/tt ~ 0.4502. 
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FIG. 3. The dependence of the Luttinger liquid quantum 
parameter A, Eq.(17) on De Boer's number Ao, Eq.(4). The 
dashed line is A = Ao. 

The A(A ) dependence found by inverting Eq.(38) is 
shown in Fig. 3 by a solid line. In the classical limit, 
A — > 0, we have A — > Ao as expected; then the effective 
pair potential (34) reduces to its bare Morse form. As 
Ao grows, the difference between A and Ao increases. We 
note that De Boer's number Ao (4) is always smaller than 
its Luttinger liquid counterpart A (17) which reflects soft- 
ening of the Luttingcr liquid by zero-point motion. The 
end point of the A(Ao) dependence which is the limit of 
stability of the liquid phase is a critical phenomenon - 
there the A(Ao) dependence has an infinite slope. Upon 
approaching Aq s from below we find that A is given by 



A = A s - 



2e 



ln(l + tt/2) 



(A 0s -A )] 1/2 



(39) 



Although both the reduced rms fluctuation /* (35) and 
bond length Q (36) are increasing functions of the Lut- 
tinger liquid parameter A (17), the relative fluctuation 
f*/Q has a maximum at Aln(l+7r/2) = Q /3. While the 
discussion of specific molecular substances will be post- 
poned until Section VI, here we note that in our attempts 
to fit real two-body potentials into the Morse form the 
reduced classical bond length was always found to satisfy 
Qo > 5. The relatively large value of Qo is a reflection 
of the strength of the short-distance overlap repulsion 
and weakness of the large-distance attraction in the pair 
interaction potential. With Q — 5 (used hereafter for 
estimates) we find that A= (5/3) ln _1 (l + tt/2) ~ 1.7651 
at the maximum of the relative fluctuation. But the Lut- 
tinger liquid cannot exist for A > X s = ln _1 (l + 7r/2) ~ 
1.0591. 
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Therefore the relative fluctuation, f */Q, reaches its 
maximal value at the border of existence of the Lut- 
tiner liquid, A s ~ 1.0591, with the magnitude not ex- 
ceeding the level of about 0.18. This fact resembling 
Lindcmann's empirical criterion of melting [34] verifies 
that our description of the bond dynamics as due to cou- 
pling to a bath of harmonic oscillators is quantitatively 
correct. The inequality f*/Q <C 1 implies that under- 
lying particles never come into close contact with each 
other - the effect of particle statistics is negligible. In 
addition, having Qo substantially larger than unity jus- 
tifies the nearest-neighbor interaction approximation in 
our treatment of the liquid. 




FIG. 4. The reduced energy per particle of the Luttinger 
liquid phase as a function of De Boer's number Ao (4). The 
lower grey scale curve corresponds to approximate solution 
described by Eqs.(26) and (27) while the upper curve is an ac- 
curate solution given by Eqs.(37) and (38). The arrow shows 
the location of the dimer dissociation threshold A02 = v2/ir. 

The only remaining approximation which needs to be 
addressed is our replacement of the vibrational spectrum 
of the system, ui{k) = (2c/H)\sm(kH/2)\, by the De- 
bye model, Lo(k) = c\k\, valid for to < lod = ttc/H. 
Both spectra are fairly close to each other: they coin- 
cide in the long- wavelength limit, kH <C 1, and end at 
the edges of the first Brillouin zone, k = ±ir/H. Physi- 
cally the chain with the Debye spectum is less susceptible 
to short- wavelength fluctuations than the original system 
- the Debye frequency lob is tt/2 times larger than the 
maximal allowed frequency. This difference necessary to 
have the correct number of degrees of freedom will only 
have a marginally small effect on final results because of 
the dominant role played in one dimension by low-energy 
long-wavelength fluctuations where the Debye approxi- 
mation becomes exact. 
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FIG. 5. The reduced quantunrtexpansion Q — Qo (a), and 
the Debye temperature 6* (b) of the Luttinger liquid as func- 
tions of De Boer's number Ao (4). Various molecular sub- 
stances are shown as solid dots on the curves. The numbers 
in the brackets on the quantum expansion graph are classical 
reduced bond lengths Qo while the size of vertical bars equals 
the relative bond fluctuation. The numbers in the brackets on 
the Debye temperature graph are Luttinger liquid exponents 
9 (11)- 

The upper curve of Fig. 4 shows the reduced ground- 
state energy per particle of the Luttinger liquid as a 
function of De Boer's number Ao (4) found by combining 
Eqs.(37) and (38). Qualitatively similar dependence was 
found in Ref. [19]; our energy is about 10 percent lower. 
For the purpose of comparison the lower grey scale curve 
of Fig. 4 shows the approximate reduced energy per par- 
ticle, Eqs.(26) and (27), which was argued to constrain 
the -E£ L (Ao) dependence from below. The quantitative 
difference between the curves is not very large and be- 
comes noticeable only for Ao > 0.15. Since variational 
analysis always constrains the ground-state energy from 
above, the true 2££ L (Ao) dependence must be sandwiched 
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between the curves of Fig. 4. We expect however that 
E*I l {\q) is well-approximated by Eqs.(37) and (38) (the 
upper curve of Fig. 4). 

Variation of other properties of the Luttinger liquid 
with De Boer's number Ao (4) can be readily calculated. 
The expression for the reduced sound velocity c* = c/c 
can be found to be 

c* - AoQ/AQo 

= e- x ln d+-/2) [! + 3(A/Q ) ln(l + tt/2)] , (40) 

which together with (38) parametrically determine 
c*(Ao). In the range of interest Qo > 5 the reduced 
sound velocity is a monotonically decreasing function of 
Ao . At the boundary of existence of the Luttinger liquid 
the reduced sound velocity reaches its minimal (and fi- 
nite) value and has a square-root singularity implied by 
Eq.(39). 

Similarly the dependence of the Luttinger liquid expo- 
nent g (11) on Ao (4) can be found by combining Eqs.(19), 
(36) and (38). Again for Qo ^ 5 the correlation exponent 
g is a monotonically increasing function of Ao reaching its 
maximal value at the boundary of existence of the liquid 
phase. 

The A dependence of the reduced Debye temperature 
9* is determined by 

9* =2K/v = 2-K 2 \e- 2X ^ 1+ «l 2 \ (41) 

combined with Eq.(38). The right-hand side of (41) has 
a maximum at A = A s /2 = 0.5 In -1 (1 + tt/2) ~ 0.5295 
which is inside the range of existence of the liquid phase. 
In view of (38) this corresponds to De Boer's number 
A = l/pe 1 / 2 ln(l + tt/2)] ~ 0.3212 which is the location 
of the maximum of the 9*{\ ) dependence. 

Fig. 5 shows the dependence of the reduced quantum 
expansion Q—Qo and Debye temperature 0* on De Boer's 
number A . The quantum expansion is found by com- 
bination of Eqs.(36) and (38); it is merely a magnified 
A(A ) curve of Fig. 3. The 9*(Xo) dependence (Fig. 5b) is 
constructed by combining Eqs.(41) and (38). Both prop- 
erties are linear functions of De Boer's number in the 
Ao — > limit. 

D. Zero-pressure phase diagram 

In order to construct the phase diagram of the sys- 
tem we need to compare all the branches of the reduced 
energy function (5); for given De Boer's number Ao the 
branch with lowest E*(Xo) singles out the ground-state. 
The outcome is shown in Fig. 6 where the reduced energy 
per particle of the Luttinger liquid, Eqs.(37) and (38), 
is drawn together with those for diatomic, Eq.(7), and 
monoatomic, E^ mo (Ao) = 0, gases. 



We see that as Ao increases away from the classical 
limit, Ao = 0, the reduced energy per particle of the Lut- 
tinger liquid increases, and at \ 0ev ~ 0.3365 a crossing 
with the dimer gas energy curve, Eq.(7), takes place - 
liquid phase evaporates into diatomic gas via discontin- 
uous transition. We note that the Luttinger liquid can 
still coexist (as a metastable state) with gas phases in 
the narrow range 0.3365 < A < 0.3896. This is shown 
by the grey scale part of the energy curve. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

0.1 I . . 1 1 . 1 1 . r- 




FIG. 6. The dependencies of the reduced energy per parti- 
cle on De Boer's number Ao (4) for various phases. The bold 
parts of the curves correspond to ground states of the system 
while the grey scale segments indicate metastable states. The 
loci of a series of substances are shown by solid dots, and the 
arrow pointing down is the dimer dissociation threshold. 

The diatomic gas is the ground-state of the system in 
the range 0.3365 < A < \/2/7r ~ 0.4502: at A 02 = V^/tt 
it undergoes a continuous dissociation transition into a 
monoatomic gas. For Ao > V2/tt quantum fluctua- 
tions are too strong and no bound state can exist: a 
monoatomic gas is the only possible state of the system. 

In determining the ground-state one needs to take 
into account all possible competing phases of the system 
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in question. We however neglected the possibility that 
trimcrs (and generally TV-atomic molecules) may come 
into play. This is because we only know how to treat the 
N = l,2,oo cases in a controlled fashion. Although it 
appears unlikely, we cannot rule out that for sufficiently 
large Ao, N > 2-atomic gases might become relevant; 
resolving this issue is left for future study. 

A critical reader still may argue that our prediction of 
the diatomic gas phase is an artifact of the variational 
treatment - the exact energy curve corresponding to the 
Luttinger liquid may go lower then what Fig. 6 shows. If 
this is the case, then the dimer gas ground state may 
disappear altogether. 

A finite-pressure treatment described next provides ad- 
ditional evidence that physics is incomplete with only 
liquid and monoatomic gas present. 

V. FINITE PRESSURE 

At zero pressure there is a qualitative difference be- 
tween a liquid which is a bound many-body state and a 
gas which is a collection of infinitely far separated par- 
ticles. Arbitrarily small confining pressure necessarily 
brings a gas to a finite density. As a result the two gas 
phases previously discussed turn into Luttinger liquids. 
The difference between the "parent" Luttinger liquid and 
what used to be a monoatomic gas becomes merely quan- 
titative - they will have differing densities, sound veloc- 
ities, correlation exponents and other properties. The 
liquid of dimers is more complicated as in addition its 
oscillation spectrum will have an extra optical branch. 

To avoid confusion we will keep referring to these 
pressure-induced Luttinger liquids as gases. As the pres- 
sure and Dc Boer's number Ao change, they may undergo 
gas-gas and gas-liquid transitions. At sufficiently large 
pressure and Ao the quantitative difference between the 
gases and the liquid must disappear. 

In what follows we will not be able to discuss these 
transitions as at the moment it is unclear how to de- 
scribe the effect of pressure on the diatomic gas phase in 
a controlled fashion. On the other hand, generalization 
of our formalism to monoatomic gas and Luttinger liq- 
uid is straightforward. Thus we will restrict ourselves to 
finding ranges of existence of these phases. 

The difference between the liquid and gas becomes 
most extreme at negative pressure. Here the liquid may 
still exist as a metastable state while the gas phase is 
impossible. 

In one dimension the pressure has dimensionality of a 
force. Assume our system is compressed by a constant 
force p applied to its ends. The system responds by ex- 
erting an outward force of magnitude p on the compress- 
ing agent which corresponds to the definition of positive 
pressure. Similarly, if the system is stretched by an exter- 
nal force, it responds by exerting an inward force on the 



stretching agent which corresponds to the definition of 
negative pressure. These two cases will be distinguished 
by the sign of p. 

Since every particle of the system is in mechanical equi- 
librium, the whole effect of pressure translates into re- 
placing the bond potential V(h) by V(h) + ph, the to- 
tal potential energy in the external field [23]. Then the 
finite-pressure analog of Eq.(20) becomes 

E*(j, Q) = v- 1 ln(l + 7) - 2e Q °- Q (l + 7" 1 )* 

+ e 2(Qo-Q) (1 + 7 -l)4A +p *g ; (42) 

where p* = pl/e is the reduced pressure. Eq.(42) should 
be minimized with respect to 7 and Q, and then v = 
X/ttXq and 7 = 2/ir substituted in the outcome will guar- 
antee translational invariance. The results are finite pres- 
sure analogs of Eqs.(34), (36), and (38): 

E*(Q) = (ttAqVA) ln(l + 2/tt) - 2e Q °- Q+xh ^ 1+7 '^ 

+ e 2(Q -Q)+4Aln(l+^/2) +p *g ) ^ 

Q = Qo + 3Aln(l+7r/2) 

- ln[(l + v / l + 2p*e 2A1 "(i+^/2)) /2] , (44) 

Ao = Ae- Aln ( 1+OT / 2 )[(1 + 2p*e 2Xln(1+ */V 

+ v /l + 2p*e 2Aln ( 1 + 7r /2))/ 2 ] 1 /2 (45) 

The expressions for the reduced rms fluctuation (35) and 
Debye temperature (41) as functions of A remain the 
same while the A(Ao) dependence is determined by (45). 
Similarly the reduced sound velocity is given by the first 
representation of Eq.(40) with Q and A determined by 
Eqs.(44) and (45), respectively. 

We note that the parameter A (17) now accounts for 
both the effects of pressure and zero-point motion. 

A. Classical limit 

In the classical limit Ao , A — ► and p* > the posi- 
tion of the minimum of (43) given by (44) naturally shifts 
to values smaller than Qo. In addition Eq.(44) predicts 
that at a very large pressure Pq ~ 2e 2Qo = 2e 10 ~ 44000 
the bond length vanishes. This conclusion is an arti- 
fact because for small interparticle separation the Morse 
potential underestimates the true strength of overlap re- 
pulsion - the bond length can only go to zero in the limit 
of infinite pressure. This flaw implies that only p* <C Pq 
results are credible which is not really restrictive as p$ is 
unrealistically large. 

The pressure dependence of the reduced sound velocity 
follows from the first representation of (40), and Eqs.(44) 
and (45). 

For not very large negative pressure the pair potential 
(43) has a minimum given by (44) and a maximum - the 
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"broken" ground state of the chain is separated from the 
metastable stretched crystal by a potential barrier. As 
the magnitude of the pressure increases, the amplitude 
of the barrier decreases, and at the classical limit of me- 
chanical stability, p* ms = —1/2, the barrier disappears al- 
together - no bound state can exist for p* < —1/2. As the 
pressure approaches the limit of mechanical stability from 
above, the reduced bond length tends to Q = Qo + In 2 
while the reduced sound velocity vanishes as (1 + 2p*y/ 4 
in agreement with general arguments of Ref. [23] . 



B. Quantum case 

With quantum effects included, the pressure at which 
the bond length vanishes (and the Morse potential ap- 
proximation fails) will be even larger than its classical 
counterpart as zero-point motion counteracts the com- 
pression. Therefore the restriction p* <C p$ remains un- 
changed in the quantum case. 

Since quantum fluctuations lead to the expansion and 
softening of the liquid phase, a smaller in magnitude neg- 
ative pressure will suffice to destabilize the liquid - the 
dependence of the limit of mechanical stability on De 
Boer's number (4), £>„ s (Ao), should be a monotonically 
increasing function of A satisfying p^ s (0) = — 1/2 (clas- 
sical limit) and p^ ns (\o s ) = (zero- pressure limit of sta- 
bility of the Luttinger liquid) . 

0.6 




FIG. 7. The A (A) dependency Eq.(45), for suffifciently 
small positive pressure (p* = 0.002 is shown). The lines of 
fixed Ao = Aoa,6 are the limits of stability of gas and liquid, 
respectively. 

For p* = the right-hand-side of (45) vanishes both 
at A = and A = oo reaching a maximum at A = X s 
which determined the limit of stability of the liquid in the 
zero-pressure case (see Section IVC). For finite positive 
pressure and A — > oo the right-hand-side of (45) behaves 
as p* 1 / 2 \ which implies that for not very large p* the 



Ao(A) dependence is a nonmonotonic function which has 
both a maximum and a minimum. The position of the 
minimum shifts to infinity as p* — > 0. An example of the 
Ao(A) dependence, Eq.(45), for sufficiently small positive 
pressure is displayed in Fig. 7 where we also show two 
lines of constant A to help identify possible phases of 
the system. 

For fixed small positive pressure and Ao < A 0a Eq.(45) 
has a unique solution for A describing the liquid. For 
Aoa < Ao < Aob Eq.(45) has three solutions. Out of 
them only the smallest (corresponding to liquid) and the 
largest (corresponding to gas) are physical. For A > A b 
there is only one solution for A describing a gas phase. 
The liquid and gas phases can coexist in the range of Ao 
between the limit of existence of the gas phase, A 0a , and 
that of liquid, Aoh- 

If the condition 



p* cxp[(2A /p* 1/2 ) ln(l + tt/2)] » 1 



(46) 



holds, then the explicit Ao-dependence of the properties 
of the gas phase can be deduced from Eqs.(44) and (45) 



A = A /p 



*l/2 



Q-Qo 
1 



(l/2)]n(p*/2) 

71" N 



IV 



-(A /p* 1/2 )in(i+i)i 



(47) 



(48) 



We note the range of applicability of these results is 
rather wide - small pressure and nonzero Ao, large pres- 
sure and arbitrary Ao, and arbitrary pressure and large 
A . The pressure dependence of the reduced length per 
particle Q (48) is the equation of state of the Morse gas. 

As the pressure increases, the distance between the 
minimum and maximum of the right-hand-side of (45) 
decreases, and at the critical pressure p* = p* ~ 0.0185 
the difference between the properties of liquid and gas 
disappears for the first time. At this pressure and Ao c — 
0.4387 the size of the liquid-gas coexistence region shrinks 
to a point. 

The negative pressure analysis is similar to what we did 
for p* = in Section IVC. For negative pressure of suf- 
ficiently small magnitude the right-hand-side of Eq.(45) 
vanishes both at A = and l + 2p*e 2Aln ( 1+7r / 2 ) = reach- 
ing a maximum in between. If Dc Boer's parameter Ao is 
below this maximum, then Eq.(45) has two solutions for 
A. The smaller (physical) solution describes a metastable 
Luttinger liquid. As A increases, the two solutions ap- 
proach each other. When Ao reaches the maximum of the 
right-hand-side of (45), we are at the limit of mechanical 
stability of the system - no liquid can exist for larger A . 

Alternatively, for sufficiently small fixed A the height 
of the maximum of (45) decreases upon increase of the 
magnitude of pressure, and at some Pm S (A ) the max- 
imum of (45) reaches the level of Aq thus bringing the 
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system to the limit of mechanical stability. It is curi- 
ous that in the quantum case the "liquid" solution dis- 
appears before the condition 1 + 2p*e 2X ln (i+*72) = q 
is reached. Therefore at the stability threshold both 
the energy barrier (between "broken" ground-state and 
stretched metastable liquid) and sound velocity remain 
finite. Only in the classical limit A — > do these quan- 
tities vanish. 

The pressure-De Boer's parameter diagram showing 
ranges of existence of liquid and gas is displayed in Fig. 8. 




FIG. 8. The p*(Xo) diagram showing ranges of existence 
of liquid and monoatomic gas phases. Different scales are 
selected on the positive and negative parts of the pressure 
axis. The liquid and gas can be in equilibrium along the 
dotted line. 

We deliberately selected different scales on the positive 
and negative parts of the pressure axis in order to be 
able to show the complete picture. As a result of this 
choice there is an illusory change of slope of the line of 
mechanical stability of liquid at zero pressure - in reality 
the J3„ s (Ao) dependence is smooth. 

The point C having coordinates Ao c — 0.4387, p* ~ 
0.0185 where the limits of existence of liquid and gas meet 
is a candidate for the liquid-gas critical point. Then the 
line of a liquid-gas evaporation transition should also pass 
through C. This curve, found by equating the ground- 
state energy (43) for both phases is shown in Fig. 8 by a 
dotted line. It does not end at C, and everywhere within 
the DC segment of the metastability line the Luttinger 
liquid has lower energy than the gas. These results imply 
that a direct liquid-monoatomic gas transition is impos- 
sible. Another gas phase, diatomic, must intervene. Al- 
though the existence of this phase will set phase bound- 
aries at finite pressure, it will not affect the ranges of 
existence of the Luttinger liquid and monoatomic gas. 




FIG. 9. The dependence of tlAgeduced quantum expansion 
Q — Qo on De Boer's parameter Ao at a negative pressure of 
large magnitude (a), within the liquid-gas coexistence region 
(b), and past the point C of Fig.8 (c). 
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Other finite pressure properties of the system can be 
readily found. As an example, Fig. 9 illustrates the depen- 
dence of the reduced quantum expansion Q — Q on Dc 
Boer's number Ao for three different pressures. The curve 
of Fig. 9a shows the expansion of the Luttinger liquid due 
to the combined effect of zero-point motion and negative 
pressure of sufficiently large magnitude (p* = —0.1 was 
used). The end point of the dependence is the limit of 
mechanical stability of the liquid. In Fig. 9b the pressure 
is selected to be within the region of the liquid-gas coex- 
istence (p* = 0.002 was used). As a result the quantum 
expansion dependence is not unique: the lower (larger 
density) curve corresponds to the liquid while the upper 
(lower density) curve curve is for the monoatomic gas. 
The end of the "liquid" curve and the beginning of the 
"gas" curve are limits of existence of these phases. In 
Fig. 9c the pressure is selected to satisfy the condition 
p* > p* c i.e. past the point C of Fig.8. Now there is 
no quantitative difference between liquid and gas. The 
presence of a relatively steep part around Ao — 0.44 is 
the effect of the proximity of p* = 0.02 to p* ~ 0.0185. 

To test the limits of applicability of our theory we 
also investigated the reduced sound velocity and relative 
bond fluctuation at various pressures. We found that 
at fixed Ao the reduced sound velocity initially increases 
with pressure but then at pressures exceeding the level 
of about 6000 or larger, it begins to decrease vanishing 
at the point where the bond length Q vanishes. The fall 
and vanishing of the sound velocity do not correspond 
to physical reality and have their origin in inadequacy of 
the Morse potential approximation at small interparti- 
cle distances. Similarly, the relative fluctuation diverges 
upon approaching the nonphysical Q = point. These 
artifacts do not pose practical limitations to our theory 
because they occur at unrealistically large pressures. We 
verified that if we limit ourselves to pressures not exceed- 
ing 5000, then the relative fluctuation is smaller than 0.5, 
and the behavior of other properties of the system is in 
agreement with physical expectations. 

TABLE I. Morse parameters for a series of molecular sub- 
stances and some of their computed properties at zero pres- 
sure. Substances are arranged in the order of decrease of their 
De Boer's number Ao. Blank entries correspond to gas ground 
states when interparticle separation is infinite. 
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159.6 


D 2 


4.028 


32.2 


.6900 


3.440 


.1994 


.7178 


5.703 


376.6 


Ne 


20.18 


35.6 


.5200 


3.110 


.1124 


.3591 


6.340 


632.6 


Ar 


39.95 


120 


.6600 


3.860 


.0343 


.1005 


5.949 


883.9 



VI. APPLICATIONS AND DISCUSSION 

Fitting two-body potentials of molecular substances 
into the Morse form can provide us with the depth of 
the potential well e, the interaction range I, and the po- 
sition of the potential minimum Hq. Supplemented by 
the masses of the underlying particles, this information 
is an input of our theory which then allows to determine 
De Boer's quantum parameter Ao (4), the ground state 
and virtually any property. 

The fitting procedure can introduce uncertainties be- 
cause real two-body interactions do not have the Morse 
form. We already know that the Morse potential under- 
estimates the strength of the overlap repulsion at short 
distances but as long as the reduced pressure does not 
exceed 5000, this flaw is practically irrelevant. 

The Morse potential also underestimates the magni- 
tude of Van der Waals attraction at large distances. We 
found that in the condensed state this shortcoming can 
be kept under the control by carrying out a Morse po- 
tential fit in a range of interparticle distances followed by 
a consistency check verifying that the segment of most 
probable particle location (formed by computed equilib- 
rium bond length plus/minus its rms fluctuation) is well 
inside the fitting range. 

The first four columns of Table I represent Morse pa- 
rameters of various molecular substances which are used 
to compute De Boer's number Ao shown in the fifth col- 
umn. The remaining three columns are reduced quantum 
expansion Q—Qo, bond length Q, and the energy per par- 
ticle E* (all at zero pressure) calculated using the theory 
developed in this paper. The loci of all these substances 
are also indicated on the reduced energy curve E*{\q) 
of Fig. 6. For the substances whose ground state is the 
Luttinger liquid, Fig. 5a also shows the magnitudes of the 
reduced classical bond length Q and its relative fluctua- 
tion f*/Q, while Fig. 5b gives the values of the Luttinger 
liquid exponent g. 

In computing these properties we also assumed that 
the three-dimensional form of the interaction does not 
change upon one-dimensional confinement of the parti- 
cles and that translational symmetry is preserved. Both 
these assumptions are approximations if the confinement 
is achieved in carbon nanotube bundles because inter- 
particle interaction is mediated by the carbon environ- 
ment [35] while the axial motion takes place in a peri- 
odic potential [15]. The former effect generally weakens 
interparticle attraction at large distances thus making 
the system more quantum. On the other hand the ex- 
ternal periodic potential due to the carbon environment 
has an opposite effect leading to upward renormalization 
of the mass. Therefore the properties of strictly one- 
dimensional matter may differ qualitatively from those 
of the matter inside nanotube bundles. The effect of an 
axial periodic potential may be even more dramatic, and 
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the effective mass approximation insufficient if the cor- 
rugation is strong enough to introduce a commensurate- 
incommensurate phase transition [36]. This potentially 
important effect is beyond the scope of our theory and 
cannot be discussed here. 

Before considering individual substances it is useful 
to look at the properties of one-dimensional matter as 
a whole and compare them with those of laboratory 
substances. The corresponding states analysis of three- 
dimensional molecular matter is based on the Lennard- 
Jones pair potential and except for the vicinity of the 
classical limit A = 0, it is empirical [2-4]. 

The substances in Table I are arranged in the order of 
decrease of their De Boer's number Ao which is naturally 
the same as in three dimensions. The main qualitative 
difference from the ordinary substances occurs because 
of the dominant role played by zero-point motion which 
in one dimension forbids the crystal ground state. Wc 
find that spin-polarized isotopes of hydrogen (hydrogen 
H j, deuterium D j, and tritium T f) and 3 He are 
monoatomic gases, 4 He is diatomic gas, while molecular 
hydrogen and heavier substances are Luttinger liquids. 
If we view the Luttinger liquid as the counterpart of the 
crystal in three dimensions, then the bold part of the 
reduced energy curve, Fig. 6, closely resembles its three- 
dimensional counterpart [3]. There is a change of slope 
somewhere between H 2 and 4 He which in our case is dis- 
sociation of the Luttinger liquid into a diatomic gas while 
in the three-dimensional world it is a melting transition. 

Our dependence of the reduced Debye temperature 9* 
on De Boer's number Ao, Fig. 5b, also looks very similar 
to its three-dimensional counterpart [2] , and even empir- 
ical values of the reduced Debye temperature are close 
to their computed one-dimensional analogs. We also find 
that #*(Ao) dependence has a maximum somewhere past 
molecular hydrogen; from empirical data it seems impos- 
sible to tell whether this effect is present or not in three 
dimensions. 

Before comparing quantum expansion in one and three 
dimensions, we note that the quantum theorem of the 
corresponding states applied to the Lcnnard-Jones sys- 
tem predicts that the reduced volume per particle (length 
in one dimension) Q is only determined by De Boer's 
quantum parameter Ao. At the same time for the Morse 
system the analogous statement is valid for the reduced 
quantum expansion Q — Qo- However the inspection of 
Table I (see also Fig. 5b) shows that the values of the re- 
duced classical bond length Qo = Hq jl belong to the rela- 
tively narrow interval roughly between 5 and 6.3. There- 
fore the variation of Qo from substance to substance can 
be ignored and within experimental error our results can 
be compared to their empirical counterparts [2,3]. Again 
we find that the Ao dependencies of the reduced quantum 
expansion in one and three dimensions are qualitatively 
similar. 

In three dimensions all these properties can be com- 



puted perturbatively in the Ao — > limit [2,3] with the 
conclusion that to leading order E* (Ao) — E*(0), Q — Q , 
and 9* all vanish linearly with Ao . This behavior is iden- 
tical to the Ao — > limit of our theory. 

As a final comment, we note that from the viewpoint 
of their electron transport properties all the molecular 
substances are normally insulators as they have com- 
pletely filled electronic shells. However at sufficiently 
large pressure when electron wave functions of neighbor- 
ing molecules overlap considerably, any substance should 
turn into a metal [37]. At that point our "molecular" 
approximation describing many-body physics in terms of 
additive two-body interactions fails. Typically this hap- 
pens at a very large pressure, and a different approach 
explicitly accounting for the dynamics of the electron de- 
grees of freedom is necessary. This complex problem is 
beyond the scope of the present paper. For the case of 
molecular hydrogen, however, we will be able to estimate 
this critical pressure when the metal-insulator transition 
takes place without leaving the framework of our theory. 

In subsequent discussion of individual substances we 
first present the results based on the pair interaction po- 
tential in free space. These conclusions are robust. On 
the other hand, our comments about the properties of 
matter inside nanotubes are speculative as they rely on 
the assumption that the effect of carbon environment can 
be accommodated within the framework of our theory 
by adjusting De Boer's parameter Ao- It is also impor- 
tant to keep in mind that inside nanotubc bundles there 
will be an additional interaction between different one- 
dimcnsional channels filled with absorbed substances. 
This interaction is responsible for exotic crossover effects 
which can be viewed as an effective change of space di- 
mensionality [38] . These effects are also beyond the scope 
of our one-dimensional theory. 

A. Spin-polarized hydrogen and its isotopes 

The pair interaction between two particles of the spin- 
polarized hydrogen family has been computed by Kolos 
and Wolnicwicz [39]. Ettcrs, Dugan and Palmer [40] have 
found a very good Morse fit to the Kolos- Wolniewicz po- 
tential; the Morse parameters shown in Table I are their 
values. 

Compared to the other elements in Table I, these sub- 
stances have the shallowest potential well which is only 
6.19K deep. Combined with its smallest mass, this makes 
H | the "quantummost" element with A = 0.9136. 
Spin-polarized deuterium D f is second in line with 
Ao = 0.6463 while spin-polarized tritium T j takes 
the third place, Ao = 0.5281. All these elements arc 
monoatomic gases at zero pressure as can be seen from 
Fig. 6. On the other hand at zero pressure in three di- 
mensions the heaviest of the family, T f, forms a liquid 
while H | and D ] are gases [40] . 
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External pressure confines these gases to a finite den- 
sity but because De Boer's numbers are larger than 
Ao c — 0.4387 corresponding to point C of Fig. 8, applying 
pressure is not going to turn them into liquids. 

B. Helium 

The pair interaction between two helium atoms is ac- 
curately described by the semi-empirical Aziz potential 
[41] which is 10.8K deep; this is the second entry in Ta- 
ble I. The authors of Rcf. [19] proposed the Morse fit of 
the Aziz potential with the parameters I = 0.5828A and 
Ho = 2. 89 A claiming that "the integrated square of the 
deviation of the fit from the Aziz potential does not ex- 
ceed 1% in the range of localization of a He atom" . These 
parameters produce Ao ~ 0.41 for A He which according 
to our theory makes it a diatomic gas. It is indeed ex- 
perimentally known [42] that 4 He can form very large 
dimers with the bond length of 52A. We verified how- 
ever that at interparticle distances that large the Morse 
fit proposed in Ref. [19] is very poor. 

Our own attempts to improve the fit increased the 
value of Ao bringing it into a narrow vicinity of the dimcr 
dissociation threshold A02 = \/2/ 7r, and without extra 
knowledge we could not make a decision whether Ao is 
larger or smaller than A 2- We resolved this dilemma 
by invoking the experimental result [42] that the binding 
energy of the 4 He dimcr is —\.\mK. Halving this value 
and dividing the outcome by 10.8i"C, the depth of the 
He — He potential well, produces the last entry in Table 
I, the reduced energy per particle of the diatomic gas. 
This can be substituted into Eq.(7) to recover the fifth 
entry, Ao = 0.4456. As expected, this is only marginally 
smaller than the dimer dissociation threshold A02 • Using 
the definition of Dc Boer's number, Eq.(4), we can now 
recover the interaction range (the third entry in Table 
I) to be I = 0.5350A. Finally the position of the mini- 
mum of the Morse potential H n = 2.980A was chosen to 
optimize the fit. 

Upon application of pressure the gas of 4 He dimers will 
turn into a Luttinger liquid; calculation of the properties 
of this dimer liquid is beyond the scope of our method. 

Dc Boer's quantum parameter for 3 He can be obtained 
from that for 4 He by invoking their mass ratio. This 
gives us the value quoted in the fifth column of the 3 He 
row in Table I. It is higher than the dimer dissociation 
threshold A02 thus ruling out earlier prediction [43] that 
3 He can form a dimer in one dimension. A many-body 
system of 3 He particles in one dimension will form a 
monoatomic gas with properties close to those of spin- 
polarized tritium. Similar to T j, the 3 He gas will not 
condense under pressure. 

Our result that A He forms a diatomic gas strictly in 
one dimension is in variance with earlier work [15,16,18] 
which predicted a liquid ground state with the binding 



energy of order a few to tens mK. This is the same 
order of magnitude as the energy per particle in the di- 
atomic gas. However Refs. [16,18] also predict a liquid- 
solid phase transition which is forbidden in one dimen- 
sion. 

In applying our results to nanotubes one has to bear 
in mind that 4 He atoms are strongly attracted to the 
interstitial channels inside nanotube bundles. The corru- 
gation felt by the individual atom is so strong that the 
effective mass enhancement is very large: m* ~ 18m [44]. 
There is also a weaker opposing effect: 28% reduction in 
the well depth of the pair interaction mediated by the 
carbon environment [35]. Combining these effects and 
assuming the interaction range does not change signifi- 
cantly, we find that De Boer's number will decrease by a 
factor of 3.6 away from its purely one-dimensional value 
thus implying a liquid ground state. Similar outcome is 
expected for 3 He. 

C. Molecular hydrogen 

The pair interaction between two hydrogen (or deu- 
terium) molecules is commonly described by the semi- 
empirical Silvera-Goldman potential [21]. Fig.l shows 
this potential together with its Morse fit; the calculated 
Morse parameters are quoted in Tabic I. 

As can be deduced from Figs. 5 and 6, and Table I the 
many-body system of H2 molecules is a Luttinger liquid 
with strongest effects of zero-point motion. It is charac- 
terized by —5.14K cohesive energy (ground-state energy 
per particle) which is an 84% reduction in magnitude 
away from the depth of the H2 — Hi potential, largest 
Debye temperature of 121K, largest Luttinger liquid ex- 
ponent g ~ 0.1, largest quantum expansion of 0.82A, and 
largest, just under 18%, relative fluctuation of the bond 
length. 

The equilibrium distance between the Hi molecules is 
4.26 A. As can be seen from Fig.l in a range around the 
equilibrium bond length significantly exceeding its rms 
fluctuation the Morse potential is a very good approxima- 
tion to the Silvera-Goldman potential. The fit worsens at 
interparticle separations exceeding 5A; there are also de- 
viations from the Silvera-Goldman potential at distances 
smaller than 3A. However at these compressions Hem- 
ley and collaborators [45] have found a softening effect 
unaccounted for by the Sivera-Goldman potential. Wc 
verified that the Morse potential shown in Fig. 1 provides 
a very good fit to the Hcmlcy-corrcctcd version of the 
Silvera-Goldman potential. 

Previous work [17] finds a liquid ground-state with the 
energy per particle to be — 4.8if and the bond length 
of 4.6A These values are close to our results. However 
we disagree with the existence of a liquid-solid transition 
found in Ref. [17] at higher density; such a transition is 
forbidden in one dimension. 
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The building blocks of one-dimensional molecular hy- 
drogen are H 2 molecules whose size of 0.75A [37] is sig- 
nificantly smaller than the computed intermolecular dis- 
tance of 4.26A Such structure can be understood quali- 
tatively from a complementary viewpoint: 

There is exactly one electron per every hydrogen atom, 
and if the protons are arranged equidistantly, then the 
valence band is half-full, and the resulting system is an 
alkali metal [1]. However Peierls [46] noticed that in one 
dimension the energy can be further lowered by displac- 
ing every second nucleus by a prescribed distance. As a 
result of the period doubling the valence band becomes 
full, and the resulting dimer chain is an insulator. We 
conclude that one-dimensional molecular hydrogen is an 
example of Peierls-distorted one-dimensional structure; 
the dimers are hydrogen molecules. 

Peierls' arguments rely on an adiabatic approximation 
which ignores zero-point motion of the nuclei; the former 
may change the answer qualitatively. Our theory which 
starts from interacting H 2 molecules shows that even for 
one-dimensional hydrogen the ground-state is a Peierls- 
distorted insulator despite strong quantum fluctuations. 

This conclusion may change upon application of pres- 
sure which brings hydrogen molecules closer to each other 
and increases relative fluctuation of the bond length. 
We argue that the distortion disappears and thus an 
insulator-metal transition takes place when all the hydro- 
gen atoms become translationally identical. This tran- 
sition is a one-dimensional version of the metallization 
transition predicted by Wigner and Huntington [37]. In 
one dimension the mechanism of the transition consists 
in "undoing" the Peierls distortion. 

Since the bond between the two hydrogen molecules 
is significantly softer than that holding the H 2 molecule 
together, we assume that upon application of pressure, 
only the former decreases. Thus the translational equiv- 
alence of all the hydrogen atoms will be achieved when 
intermolecular spacing reaches the value of order 1.5A, 
twice the size of the H 2 molecule. With zero-pressure 
intermolecular spacing being 4.26A, this corresponds to 
compression by a factor of 2.84. The corresponding re- 
duced pressure p* » 1 can be found by inverting Eq.(48): 

p* ~ 2e 2{Qa - Q+[V ^ Xoln{1+ i ) -i ]eQ ' Qo} (49) 

Substituting here Q = 4.986, Q = H/l = 2.174, and 
A = 0.2819 we arrive at the reduced pressure of 545. Our 
theory which does not explicitly consider electronic de- 
grees of freedom fails in the vicinity of the inverse Peierls 
transition. 

For molecular hydrogen in one dimension the unit of 
pressure is a force of e/l ~ 6.44 * lO -12 ^ strong. Mul- 
tiplying this by 545 we find that 3.51 * 1CP 9 7V force 
compressing one-dimensional hydrogen may suffice to in- 
duce a transition into metallic state. If this force is 
applied at the 1A 2 area, then the corresponding three- 



dimensional pressure will be 351GPa. We note that 
three-dimensional solid hydrogen subject to pressure that 
big still resists metallization [47]. Unfortunately the ac- 
curacy of our estimate is not great because of the expo- 
nential dependencies in (49) - the actual one-dimensional 
transition may happen at lower or larger pressures. 

For hydrogen confined inside interstitial channels of 
carbon nanotube bundles the carbon environment effec- 
tively reduces the well depth of the pair interaction by 
54% [35]. This effect alone would suffice to turn the 
many-body system of hydrogen molecules into a gas of 
(#2)2 complexes. However if the effective mass enhance- 
ment is comparable to that for He [44] , the liquid ground 
state might be restored. 

As can be seen from Fig. 6 and Table I, the ground 
state of molecular deuterium is a Luttingcr liquid. The 
quantitative difference from the properties of molecular 
hydrogen is solely due to the fact that D 2 has a larger 
mass. The cohesive energy of the one-dimensional D 2 
liquid is —Yl.XK (62% reduction in the magnitude of 
the D 2 — D 2 pair potential well), the Debye tempera- 
ture is 99. 8K while the equilibrium distance between D 2 
molecules is 3.94A 

D. Heavier substances 

While discussing the physics of molecular hydrogen in 
one dimension we came to the conclusion that it can be 
viewed as an example of Peierls-distorted structure. The 
same arguments are applicable to any element with odd 
number of electrons: the period doubling should take 
place and the resulting system must be an insulator. We 
note that three-dimensional counterparts of these sub- 
stances are metals. These observations imply that our 
theory is also applicable to substances which tradition- 
ally are not considered to belong to the molecular group. 
For example, one-dimensional lithium must be an insu- 
lating Luttinger liquid of Li 2 molecules [48]. Similar to 
molecular hydrogen, under pressure it should undergo a 
metallization transition. Had we known the pair interac- 
tion between two Li 2 molecules, we could have computed 
the properties of the lithium liquid. From the viewpoint 
of the quantum theorem of corresponding states molec- 
ular lithium is expected to occupy a place somewhere 
between D 2 and Ne. 

Table I also contains the Morse data for Ne and Ar. 
They were obtained from the parameters of the Lennard- 
Jones interaction potential [1] which is commonly used 
to describe these noble gases. Some of the properties 
of these substances in one dimension can be found in 
Table I and in Figs. 5 and 6. The equilibrium interparticle 
spacing, cohesive energy and Debye temperature can be 
extracted from what is shown in the same manner as was 
done for lighter elements. As the underlying particles 
become heavier, the effect of zero-point motion decreases. 
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For elements heavier than Ar quantum fluctuations can 
be ignored for most practical purposes. 
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